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ABSTRACT 

We report the detection of the "''^Sc nuclear decay hnes at 67.9 keV and 78.4 keV associated with the 
nuclear decay of ^4^^ y^ Cassiopeia A. The line emission was observed by the PDS instrument on board 
BeppoSAX, which recently observed the supernova remnant for over 500 ks. The detection of the line 
emission with a flux of (2.1 ±0.7) 10~^ ph/cm^/s in each line (90% confidence) is at the ba significance 
level, if we can assume that the 12-300 keV continuum is adequately represented by a single power law. 
However, as the nature of the continuum is not clear we investigate various other possibilities. A more 
conservative estimate of the line flux is made by assuming that a power law continuum is at least a good 
approximation to the continuum emission for a narrower 30-100 keV energy range. With this limitation 
the measured line flux is (1.9 ± 0.9) 10~^ ph/cm^/s, with the detection still at the 3.4(t significance 
level. We suggest that together with the CGRO-COMPTEL measurement of the '''*Ca fine at 1157 keV 
of (3.3 ± 0.6) 10"^ ph/cmVs a flux for aU three lines of (2.5 ± 1.0) 10"^ ph/cmVs for Cas A can be 
adopted. This imphes an initial "''Ti mass of (0.8 - 2.5) lO^"* Mq. 

Subject headings: gamma rays: observations - ISM: individual (Cassiopeia A) - nuclear reactions, 
nucleosynthesis, abundances - supernova remnants - X-rays: ISM 



2. OBSERVATIONS AND DATA 

In May and June this year BeppoSAX observed Cas A 
for more than 500 ks. The primary goal of this obser- 
vation was to detect the 67.9 keV and 78.4 keV nuclear 
decay lines with the hard X-ray Phoswich Detection Sys- 
tem, PDS (Frontera et al. 1997). The PDS consists of four 
NaI(Tl)/CsI(Na) scintillation detectors behind two rock- 
ing collimators with a 1.3° field of view. During the obser- 
vation the collimators switch back and forth between the 
on source position and two opposite off source positions, 
such that always one collimator observes the source and 
the other the background. Due to the alternating posi- 
tion of the collimators the effective exposure of the source 
with the PDS is about half the total observation time. In 
this case the effective PDS exposure is 256 ks. We used 
additional archival PDS spectra of on source Cas A ob- 
servations, resulting in a total effective PDS exposure^ of 
311 ks. 

As explained in the BeppoSAX analysis cookbook (Fiore 
et al. 1999), the standard software pipeline incorporates 
two methods for rejection of background events due to 
energetic particles. The more advanced method, which 
is preferred for a background dominated source like Cas 
A, uses a variable, energy-dependent, rise-time rejection 
threshold. This results in a better particle rejection, but 
at the cost of a 7% loss in effective area. Although this is 
the preferred method, to get a feel of the systematic errors 
involved, we also include some analysis of the same obser- 
vations, but with spectra extracted using the constant rise- 

^ Chandra fellow 

2 Columbia Astrophysics Laboratory, Columbia University, MC 5247, 550 W 120th street, New York, NY 10027, USA 

3 Naval Research Laboratory, Code 7674L, Washington DC 20375, USA 

* Space Research Organization Netherlands (SRON), Sorbonnelaan 2, NL-3584 CA, Utrecht, The Netherlands 

^ The observational identification numbers for the additional archival data used for our analysis are 30011001, 30011002 and 30795005. They 
were chosen on the availability of spectra made with the variable rise time rejection method. 



1. INTRODUCTION 

The detection of gamma ray line emission at 1157 keV 
from the supernova remnant Cassiopeia A by CGRO- 
COMPTEL (lyudin et al. 1994) led to renewed interest 
in the nucleosynthesis and properties of '*'*Ti, the radio- 
active element with which this '^^Ca nuclear de-excitation 
emission is associated. 

However, the decay chain ^'^Ti -^ "^^Sc -^ "'''Ca produces 
two other nuclear de-excitation lines of '^''Sc at 67.9 keV 
and 78.4 keV with a flux equal to that of the 1157 keV 
emission (see Diehl & Timmes (1998) for a review). As ob- 
servations with the hard X-ray instruments CGRO-OSSE, 
RXTE-HEXTE and BeppoSAX-PDS failed to detect those 
hnes (The et al. 1996; Rothschild et al. 1997; Vink et al. 
2000), some doubt was cast on the observed 1157 keV line 
flux, if not on the detection itself. With subsequent obser- 
vations by COMPTEL, the flux estimates changed from 
(7.0 ± 1.7) 10"^ ph/cmVs to (3.3 ± 0.6) 10"^ ph/cm^/s 
(lyudin 1997), whereas a 83 ks BeppoSAX-PDS spectrum 
narrowed the ^'^Sc fluxes to a 99.7% upper limit of 4.1 10~^ 
ph/cmVs (Vink et al. 2000). 

In this paper we report on a new, deep, exposure 
of Cas A with the BeppoSAX X-ray observatory. The 
BeppoSAX-PDS spectrum from this observation, together 
with archival data, finally gives a reliable estimate of the 
the 67.9 keV and 78.4 keV fine fiuxes. 
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Fig. 1. — BeppoSAX-PDS spectrum of Cassiopeia A. The solid line in the left figure shows the spectral model including the ''''Sc lines at 
68.9 keV and 78.4 keV. The dotted line shows the expected systematic error in the background subtraction, as estimated from blank field 
PDS spectra (Guainazzi & Matteuzzi 1997). The figure on the right shows the best fit power law model without ''''Sc lines. The channels 
above 80 koV have been rebinned for presentational purposes. 



time rejection method. Compared with the statistical un- 
certainties, the measured "^^Sc flux does not depend much 
on the method chosen, but statistical errors are slightly 
larger for spectra that were extracted using the constant 
risetime method (Table 1). 

There exists a well known error in the absolute flux den- 
sity scale of the PDS instrument of a factor 0.86±0.03. We 
took this factor and the 7% effective area loss into account 
by multiplying the response matrix by a factor 0.80. Al- 
though the statistical error of the flux calibration is small, 
there exists some scatter in this relation (Fiore et al. 1999). 
In view of this, and the small differences between measure- 
ments made using the two rise time rejection methods, we 
adopt a conservative 10% systematic error. In order to 
obtain reliable x^ statistics, the spectra were rebinned to 
channel widths of approximately Ai?/3, where AE is the 
full width at half maximum (FHWM) energy resolution, 
which is about 9 keV at 75 keV. Further rebinning was not 
necessary as the data are background dominated, result- 
ing in gaussian error distributions for each channel. We 
verified that the spectra of each collimator were consistent 
with the best fit model to the overall spectrum (see be- 
low). Similarly, we verified the consistency of the recently 
observed spectra and the archival spectra. 

As far as we know there are no instrumental background 
lines coinciding in energy with the ^^Sc lines. Possible 
contamination sources of line emission close to the ^"^Sc 
lines are tantalum from the collimator and the on board 
^''^ Am calibration source. ^''^Am emits photons at 60 keV, 
whereas the Ta Ka and K/3 transitions are at 57 keV 
and 65 keV with an emission ratio of 4:1. These lines 
may be responsible for a small emission excess around 
60 keV (Fig. 1). However, adding these lines to our mod- 
els changed the measured '*'*Sc line flux with less than 5%; 
much smaller than the statistical uncertainty in the line 
flux. 

Note that the observed power law normalization of the 
PDS spectrum agrees with that of the RXTE-HEXTE 
(Allen et al. 1997). At 15.9 keV their model fits imply 
2.6 10^^ ph/s/kcV, whereas the PDS spectrum indicates 



2.510 '* ph/s/keV. The photon index above 16 keV of the 
RXTE-HEXTE spectrum, 3.04 ± 0.15, agrees within 1.6cr 
with the value reported below. The continuum normaliza- 
tion also agrees within 1.5cr with measurements by CGRO- 
OSSE (The et al. 1996), which measured a fiux density of 
(9.0 ± 2.1) lO"'^ ph/cmVkeV at 100 keV, compared to 
(5.8 ± 0.8) 10^^ ph/cm^/keV for our measurement, ob- 
tained from extrapolating the power law normalization to 
100 keV (first row in Table 1). Note that the CGRO- 
OSSE had a much larger field of view (11° x 4°4) than the 
BcppoSAX-PDS. 

3. DATA ANALYSIS 

We present the combined PDS-spectrum and the sim- 
plest possible model, consisting of a power law continuum 
and two nuclear decay lines, in Figure 1. For this mod- 
eling of the continuum shape the ^^Sc lines are detected 
with a significance of more than 5(7 (Ax^ = 30). The 
flux in each line is (2.1 ± 0.7) 10~^ ph/cm^/s and the 
photon index is 3.3 (see Table 1 for details). However, 
the nature of the hard X-ray continuum is still under de- 
bate. It has been argued that it is synchrotron radiation 
from shock accelerated electrons, or, alternatively, that it 
is bremsstrahlung emission caused by a non-thermal tail 
to the thermal electron distribution (Asvarov et al. 1990; 
Allen et al. 1997; Favata et al. 1997; Bleeker et al. 2001). 
A more specific bremsstrahlung model was worked out by 
Laming (2001a, b), who calculated the spectrum of elec- 
trons accelerated by lower hybrid waves associated with 
shocks in the ejecta. 

In this letter we limit ourselves to the detection of the 
'*'*Sc lines. We therefore leave a detailed discussion of the 
hard X-ray continuum, which should also include data for 
the energy range 0.5-10 keV, to a future article. However, 
in order to test the dependence of the measured line flux 
estimates on the assumed continuum shape we tested var- 
ious continuum models, the results of which can be found 
in Table 1. It is surprising that the best fit to the PDS- 
spectrum in a statistical sense is provided by a power law 
model plus *'*Sc line emission. However, as the measured 
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line flux depends on the continuum model chosen, we have 
to take the uncertainty about the nature of the continuum 
into account. 

It is clear from spectral fits to the narrower 30-100 keV 
spectral range that the inclusion of a thermal component 
with kTc= 4.2 keV has little effect on the estimated ^''Sc 
flux. Certainly, for this narrow range a power law is a 
reasonable approximation for the continuum. As the en- 
ergy range is smaller, there is more statistical uncertainty 
about the photon index, and therefore the ^''Sc line flux is 
more uncertain. Indeed the significance of the line emis- 
sion drops from the 5a to the 3.4cr level. The photon 
index versus line flux confidence contours for this energy 
range is shown in Figure 2. The 3a upper limit on the 
flux in both lines is 3.5 10^^ ph/cm^/s, based on the 30- 
100 keV energy range. This upper limit is comparable 
to the ^■^Ca line flux at 1157 keV recently obtained from 
CGRO-COMPTEL data, (3.3 ±0.6) 10"^ ph/cm^/s (68% 
confidence range), see lyudin (1997). However, both mea- 
surements are consistent given the 25% systematic uncer- 
tainties for the COMPTEL measurement (Dupraz et al. 
1997). 

Although it is difficult to combine the COMPTEL and 
the PDS results, due to the unknown nature of the sys- 
tematic errors, the fact that there are now two indepen- 
dent measurements of line emission associated with the 
'''^Ti decay increases the credibility of the detections. We 
therefore suggest adopting a ^^Sc/^^Ca line flux for Cas 
A of (2.5 ± 1.0) 10~^ ph/cm^/s, which is consistent with 
both the PDS and COMPTEL measurements. 

4. CONCLUSIONS 

Since the initial discovery of the '*'*Ca nuclear decay lines 
by CGRO-COMPTEL (lyudin et al. 1994), our knowledge 
of the formation and decay of ^^^j j^^^g substantially im- 
proved. Recently, the decay time of '*'*Ti has been ac- 
curately measured to be 85.4 ± 0.9 yr by three indepen- 
dent experiments (Ahmad et al. 1998; Gorres et al. 1998; 
Norman et al. 1998). As the Cas A supernova was prob- 
ably observed by the English astronomer J. Flamsteed 
in A.D. 1680 (Ashworth 1980), the age of Cas A is ~ 
320 yr. A good alternative is to use the kinematic age 
of the fast moving optical knots, 330 yr (Thorstensen et 
al. 2001), but the age difference is so small that it has lit- 
tle effect on the inferred range of initial ''^Ti mass. The 
line flux of (2.5 ± 1.0) 10^'"' ph/cm^/s, combined with the 
^''Ti decay time and a distance to Cas A of 3. 4^°'^ kpc 
(Reed et al. 1995), yields an initial ^''Ti mass in the range 



(0.8 - 2.5) 10-4 M0, with 1.2 lO^* Mq corresponding to 
the adopted ''^Sc/^'^Ca flux and distance to Cas A. 

This is rather high compared to model predictions, 
which usually indicate ^''Ti masses below 10^^ Mq, except 
for progenitor masses around 12Mo and above ~25M0 
(Timmes et al. 1996). For that reason Mochizuki et al. 
(1999) made the interesting suggestion that the "^^Ti de- 
cay may be delayed, as a result of complete ionization of 
^''Ti, inhibiting the decay of '*'*Ti by the capture of a K- 
shell electron. However, recently Laming (2001c) showed 
that the electron temperature in the ejecta was probably 
never high enough to seriously affect the ^'^Ti decay. 

The production of '^■^Ti increases with the size of the 
helium core of the progenitor, but material falling back 
on the neutron star or black hole limits the amount of 
ejected ^''Ti. Model calculations show that massive stars 
have more fall back, but pre-supernova wind loss may limit 
the amount of material falling back (Timmes et al. 1996). 
This agrees with the idea that the progenitor of Cas A was 
a not too massive (~ 3OM0) Wolf-Rayet star that suffered 
heavy mass loss (Vink et al. 1996). Additionally, assyme- 
tries in the explosion may have increased the amount of 
^''Ti synthesized (Nagataki et al. 1998). The likely pres- 
ence of a neutron star (Chakrabarty et al. 2001) in Cas A 
further constrains the explosion scenario for Cas A, as too 
much fall back would have resulted in the formation of a 
black hole. 

Although we have now finally detected the ^''Sc nuclear 
decay lines at 67.9 keV and 78.4 keV, interesting obser- 
vations remain to be done with future hard X-ray and 
Gamma-ray missions. The solid state detectors on board 
Integral will be able to measure the line broadening of the 
^''Ca line accurately, and constrain the properties of the 
hard X-ray continuum further. In addition, hard X-ray 
experiments using multi-layer mirrors will be able to map 
the spatial distribution of '^^T[ in Cas A on the arcminute 
scale. 
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Table 1 
Summary of spectral model fits. 



RT 


model 


spectral 


**Se flux 


PL 


PL norm 




radio 


roll ofl:' 


Emission 


xV- 


method 




range 
(keV) 


(10^^ ph/cm^/s) 


index 


@ 1 keV 
(ph/s/keV) 




index 


energy 
(keV) 


measure^ 

(10^^ em"^) 




V 


PL 


12 - 300 


2.1 ±0.7 


3.30 ±0.05 


2.3 ±0.3 










69.5/63 


V 


PL+therm. 


12 - 300 


1.0 ±0.7 


2.71 ±0.06 


0.28 ±0.05 








9.25'=(fixed) 


96.7/63 


V 


PL+therm. 


12 - 300 


2.0 ±0.7 


3.2 ±0.2 


1.8 ±0.8 








1.3(< 4.5) 


69.4/62 


V 


sresc-l-therm. 


12 - 300 


3.2 ±0.8 






0.85 


(> 0.848)'= 


9.2 ±0.3 


9.25''(fixed) 


120.8/63 


V 


sresc+therm. 


12 - 300 


2.8 ±0.9 






0.85 


(> 0.838)'= 


10 ± 1 


5.8 ±0.9 


75.9/62 


V 


PL 


30 - 100 


1.9 ±0.9 


3.1 ±0.4 


n-sto.? 










15.6/20 


V 


PL+therm. 


30 - 100 


1.8 ±0.9 


3.0 ±0.04 


7+^-* 

"■'-0.6 








9.25 (fixed) 


15.5/20 


C 


PL 


12 - 300 


2.5 ±0.8 


3.33 ± 0.06 


2.5 ±0.4 










74.5/63 


C 


PL 


30 - 100 


1.7± 1.0 


2.8 ±0.5 


0-30tl,il 










18.7/20 



Note. — All spectral models can be found in X-ray spectral analysis package xspec vll (Arnaud 1996). The first column indicates which 
background rejection critcrium was used, - variable (V), or constant (C) - to obtain the source spectrum (see text). The Sc lines were modelled by 
two deltalines with energies fixed at 67.9 keV and 78.4 keV, and equal line fiux. PL is short for power law. We used vmekal (a coUisional equilibrium 
model) for the thermal component (Mewe et al. 1995). The sresc model is described in Reynolds & Keohane (1999); it gives the synchrotron emission 
from a relativistic power law distribution of particles with an exponential cut-ofT. Errors are statistical errors, and correspond to 90% confidence 
intervals. 

^The emission measure is defined as —^ — ^^ . Note that in order to keep in line with models used by Vink et al. (1996), Favata et al. (1997) and 

Vink et al. (2000) we used a plasma enhanced in helium and nitrogen by a factor 10 (based on optical observations of the shocked circumstellar 
medium, e.g. Chevalier & Kirshner (1979)). The electron temperature was fixed to kT^— 4.2 keV (Vink et al. 2000). Note that the ejecta is most 
likely the main source of thermal X-ray emission in Cas A, and so an O-rich composition could also be a reasonable model. This would not, however, 
make any difference to the shape of the thermal continuum at photon energies relevant here. 

The fixed emission measure is identical to that used by Vink et al. (2000), comparing with ASCA-SISO and BcppoSAX-MECS data shows that 
this implies that below 10 keV about half of the continuum is thermal, whereas the other half comes from an additional component. 

*^ We allowed the radio spectral index, nominally 0.77 for Cas A, to vary in the range 0.7 to 0.85; the normalization parameter is the radio flux 
density at 1 GHz, which was fixed to 2720 Jy (Green 2000). The parameters derived from this model should be treated with caution, as the values 
for radio index and roll off energy are highly correlated. 



